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We examine vortex pinning and dynamics in tliin-film superconductors interacting with square and 
rectangular pinning arrays for varied vortex densities including densities significantly larger than 
the pinning density. For both square and rectangular pinning arrays, the critical depinning force 
shows maxima at only certain integer matching fields where the vortices can form highly ordered 
arrays. For rectangular arrays the depinning force and commensurability effects are anisotropic with 
a much lower depinning threshold for vortex motion in the easy fiow directions. We find evidence 
for a crossover in pinning behavior in rectangular pinning arrays as the field is increased. We also 
show analytically, and confirm with simulations, that for B — 2B^ the strongest pinning can be 
achieved for rectangular pinning arrangements rather than square for one direction of driving force. 
Under an applied driving force we find a remarkable variety of distinct complex fiow phases in both 
square and rectangular arrays. These fiow phases include stable sinusoidal and intricate pinched 
patterns where vortices from different channels do not mix. As a function of the driving force certain 
fiow states become unstable and transitions between different phases are observed which coincide 
with changes in the net vortex velocities. In the rectangular arrays the types of flow depend on the 
direction of drive. We also show that two general types of plastic flow occur: stable flows, where 
vortices always fiow along the same paths, and unstable or chaotic fiows. 



I. INTRODUCTION 

Vortex pinning in systems where defects are arranged 
in periodic arrays are an ideal system in which to study 
optimal pinning arrangements and commensurability ef- 
fects since the properties of the pinning sites, such as the 
size, periodicity of the array, and array geometry, can be 
highly controlled. These arrays can be fabricated using 
nano-lithography techniques in which periodic arrange- 
ments of micro-holes [p]-Jl2| or magnetic dots |13-19| of 
various geometries can be constructed. A more recent 
technique is the use of Bitter decoration where the mag- 
netic particles used for the initial decoration can act as 
well ordered arrays of pinning . In these systems com- 
mensurability effects can be observed in the form of peaks 
in the critical current at fields where the number of vor- 
tices equals an integer or a fraction of the number of 
pinning sites. At these commensurate fields the vortices 
can form a highly ordered lattice where vortex-vortex 
interactions which reduce the effective pinning are mini- 
mized, while at incommensurate fields the vortex lattice 
can be disordered and vortex-vortex interactions make 
the pinning less effective. Direct imaging of vortex struc- 
tures with square pinning arrays have been conducted 
with Lorentz microscopy and scanning Hall-probes 
pO| which have confirmed that the vortices form highly 
ordered crystals at the integer and some fractional match- 
ing fields. These images along with simulations [E3| have 



also shown that at the different matching fields different 
types of ordered vortex crystals can be stabilized. In the 
samples imaged by Harada et at, it was found that 
beyond the first matching field all the pinning sites are 
occupied with a single vortex and that the additional vor- 
tices sit in the interstitial regions. The particular vortex 
lattice symmetry found at the second matching field was 
square; at the third matching field there was an ordered 
lattice with alternating pairs of interstitial vortices; while 
at the fourth matching field the overall lattice was trian- 
gular. The same types of ordered vortex lattice crystals 
were also seen in simulations with square pinning arrays 
and additional types of stable vortex crystals were ob- 
served for triangular pinning arrays |p3|| . Recent scan- 
ning Hall-probe experiments [p^ have found various or- 
dered rational fractional fillings less than the first match- 
ing field and have also shown that for fields greater than 
the first matching field, vortices can have both multiple 
occupancy per pinning site as well as vortices that are 
located in the interstitial regions. 

Most experimental studies of vortex matter interact- 
ing with periodic pinning have considered square and 
triangular pinning lattices. Some recent studies, how- 
ever, have been performed with Kagome pinning arrays, 
which produced pronounced matching effects at non- 
integer matching fields . Periodic arrangements with 
rectangular geometries of magnetic dots [|l8|,|9| and holes 
| p7| , pT 22 1 have also been studied recently. In these sys- 
tems there can be two periodicities associated with the 
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two sides of the rectangular cell, a and b. In experiments 
with magnetic dots [^,|9| a field dependent crossover 
in the commensurability effects was observed. Sharp low 
field matching effects occurring at every integer matching 
field were found to cross over to much broader match- 
ing effects occurring at fields where the vortex density 
matches with the periodicity of the short side of the 
rectangular pinning array. In these systems one would 
expect an anisotropic response, with a lower pinning or 
easy flow direction through the wide end of the rectan- 
gular cell since the flow of interstitial vortices will be less 
impeded by the vortices pinned by the dots. Addition- 
ally one would still expect an anisotropic response due to 
increased vortex-vortex interactions along the short end 
of the cell at the incommensurate fields which will lower 
the pinning force for fields less than the first matching 
field, as well as for fields greater than the first matching 
field when multiple occupancy per pinning site occurs. 
Recent magneto-optical imaging experiments in samples 
with rectangular pinning arrays have found evidence for 
anisotropic vortex flow . Anisotropic pinning has also 
been observed in related systems where the individual 
pinning sites have an anisotropic geometry pO|] . 

In square pinning arrays recent simulations have shown 
that a remarkable number of distinct dynamical phases 
are possible when the number of vortices is larger than 
the number of pinning sites. These phases include one- 
dimensional flow of interstitial vortices, random plastic 
flow, soliton-likc flow along the pinning sites, and coher- 
ently moving clastic flow phases j2^]. Other simulations 
of vortices in systems with periodic pinning arrays have 
also observed transitions from different types of plastic 
flow to elastic flow p6| , p7[ | 

Experiments and simulations have also examined the 
pin ning and dynamics of vortices through thin channels 
iPSHSOj , where the vortices moving through the channels 
experience a periodic potential created by the ordered 
lattice of vortices that are immobile outside of the chan- 
nels. Peaks in the critical current are observed where 
the vortices in the channel can form a commensurate or- 
dering. Additionally simulations ||3C| ] have shown that a 
wide variety of dynamical phases occur such as soliton 
like flow, and a periodic winding motion of vortices. 

Experimentally, evidence for the flow of interstitial vor- 
tices in samples with periodic pinning arrays where in- 
dividual pinning sites are small has been observed in 
transport measurements and current-voltage curves 
Shapiro steps for the ID flow of interstitial vortices be- 
tween pinning sites with an AC and DC applied driving 
has been observed in experiments [Tl| ] and simulations 
psf . Direct imaging of vortices using Lorentz microscopy 
has demonstrated the ID flow of interstitial vortices be- 
tween pinning sites as well as pulsc-likc motion of vortices 
along pinning sites and between the pinning sites. The 
pulse or soliton flow of vortices along the symmetry direc- 
tions of the pinning arrays has also been observed with 



magneto-optical imaging . 

In systems with random pinning under an applied driv- 
ing force the vortex lattice can exhibit different types of 
transport behavior such as plastic flow, where the vortex 
lattice is highly disordered or liquid like and a portion 
of the vortices can remain immobile while other portions 
can tear past. There can also be elastic flow where the 
moving vortices keep the same neighbors and the over- 
all lattice structure can have a crystalline or smectic or- 
der. As a function of applied driving force transitions or 
crossovers can occur between different flow phases. Evi- 
dence for such transitions has been observed in transport 
measurements , voltage noise , neutron scattering 
lH, Bitter-decoration @,||], and STM ^ and sim- 
ulations [^7|-^. These moving phases have also been 
studied theoretically |^ . 

In this work we present the results of simulations for 
vortices in 2D superconductors interacting with square 
and rectangular pinning arrays. In particular we examine 
the pinning characteristics and the dynamical states. In 
previous simulations only square and triangular pin- 
ning arrays were examined and the dynamical flow states 
were only examined in detail for filling fractions less than 
2.25. In addition previous simulations calculated the 
magnetization curves but did not calculate the depin- 
ning force directly from current-voltage curves. Here, 
we examine anisotropic pinning by considering rectangu- 
lar pinning arrays for a/b = 2, 1.6, and I under applied 
driving. We have also calculated analytically the critical 
currents in the a and b direction at the second match- 
ing fields as a function of a/b and find that the critical 
current is optimal not for a/b = I but for a rectangular 
geometry when driven in certain directions. For square 
pinning arrays we find much more pronounced match- 
ing effects at integer and fractional matching fields for 
low applied fields; however, for higher fields the match- 
ing effects are reduced and certain peaks are absent, in 
agreement with recent experiments | ^ . For rectangular 
arrays with a/b — 2 we find that the shape and strength 
of the commensurability effects at the matching fields de- 
pends on whether the driving is in the a or b direction. 
The depinning force is much higher along the long direc- 
tion. We also find evidence for a crossover in the com- 
mensurability effects at high fields. Images of the vortices 
in the rectangular pinning array show that highly ordered 
vortex crystals are stabilized at the matching fields where 
a peak in the depinning force is observed, while more dis- 
ordered crystals are formed at the matching fields where 
peaks are not observed. 

Under an applied driving force we show that a re- 
markable variety of complex dynamical phases emerge 
in square and rectangular arrays. These phases depend 
strongly on the geometry of the pinning lattice. De- 
spite the complexity of these dynamic phases two gen- 
eral classes of fiow can be identified. The first is elastic 
flow of a mobile sub-lattice of interstitial vortices between 
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pinned vortices with the vortex motion in stable well de- 
fined patterns. The other type of flow is a chaotic or 
mixing flow where the interstitial vortex motion is disor- 
derly. The moving channels show a mixing effect, in that 
the vortices from one channel move to other channels. 
The stable, non-mixing flow phases occur only for certain 
integer matching fields and rational machining fields. As 
a function of applied driving force dynamical transitions 
between different types of flow states are possible which 
coincide with features in the driving force versus velocity 
curves or voltage-current curves which can be detected 
exp eriment ally. 



II. SIMULATION 

We simulate a thin superconductor where the vortices 
can be considered as 2D objects interacting with a loga- 
rithmic potential, Uy = — ln(r), where energy is normal- 
ized to Ay = ^q/SttA, $0 being the flux quantum and 
A the effective 2D penetration depth for a thin film su- 
perconductor. The normalized overdamped equation of 
motion for a vortex i is 



dt 



(1) 



The normalized force on vortex i from the other vortices 
is f™ = — X^j^i ^iUv{rij), where Uy{rij) is the effective, 
re-summed, vortex- vortex potential between two vortices 
in a computational cell with periodic boundary condi- 
tions Pinning is modeled as attractive parabolic 
wells with 
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Q is the Heaviside step function, r)f is the location 
fp is the maximum pinning force, 



(2) 



of pinning site fc. 



and Vp is the radius of the 



pinning sites. The pinning sites are placed in a rectangu- 
lar array {nL^^mLy), where n and m are integers. The 
characteristic length, by which all lengths are normal- 
ized, is chosen to be rg = ■^fp. In a typical experiment 
= 400—900 nm and the pinning radius is 125 nm. The 
characteristic time in the normalized equation of motion 
is T = rjrp/Ay, where rj is the Bardeen-Stephen friction. 
The equation of motion is numerically integrated by an 
Euler method using a normalized time step of dt — 0.02. 

The initial vortex positions arc obtained from anneal- 
ing by a high temperature, where we add a Langevin 
noise term to equation (1) such that (f/") = and 
(ff (i) ■ fj(i')) = 4T5y(5(i - t'), where T is the tempera- 
ture normalized to Ay/k^, being the Boltzmann con- 
stant. We choose an initial temperature high enough such 
that the vortices are in a molten state and cool to T = 
in 20 increments where we allow 20000 time steps be- 
tween each increment. After annealing, the driving force 



is slowly increased from /d = and the vortex positions 
and velocities are monitored using 14 = X^f'" ^ ' ^^'^ 
Vy = y ' '^i- this work we consider simulations 
where the driving is only in the x-direction or y-direction. 
The depinning force is defined as the force at which the 
vortex velocities reach 0.03/^. 



III. CRITICAL FORCE AT B = 2B, 



We first consider analytically the depinning threshold 
in the case of the vortex lattice driven in the x or y direc- 
tions at i? = 2B^ for arbitrary and Ly . At this filling 
fraction the vortex lattice consists of the vortices at the 
pinning sites along with well ordered interstitial vortices 
located between the vortices at the pinning sites. We 
make the assumption that the unpinned vortices form 
a perfect rectangular lattice, so that they effectively do 
not interact due to symmetry. With this assumption we 
need to consider only the depinning of a single intersti- 
tial vortex moving in the periodic potential created by 
the vortices at the pinning sites. The vortex will move 
one-dimensionally along the direction of the drive and 
experience the periodic potential created by the vortices 
at the pinning sites. 
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The vortex- vortex interaction term f^'"{x^y) can be cal- 
culated using the re-summation method for logarithmi- 
cally interacting particles in 2D |^^; the resulting equa- 
tion of motion for an interstitial vortex (i) moving at 
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dxi 
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^Ly along the x-direction is: 
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The maximum depinning force will occur when the 
summation term is maximum. In the limit of Ly L^, 
this happens for Xi = jL^, and the resulting critical force 



fx = ^sech(7r-^) . 



(5) 



This expression is accurate for L^^l.bLy (the error is 
about 0.2% for = Ly and about 5% for L^ = 2Ly). By 
symmetry, we trivially write the corresponding critical 
force in the y-direction, for L^ ^ Ly, to be given by. 
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which is then accurate for Ly^LbL^- 
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In figure 1 we plot the critical force in the i-direction 
and y-direction from equations (5) and (6) where Ly is 
kept at a fixed value of 2.0 and is varied from 4.0 
to 1.25. For all increasing Ly/L^ the critical force in 
the x-direction is reduced. This can be understood by 
considering that the pinning of the interstitial vortices is 
caused by the interactions with the vortices at the pin- 
ning sites. With no driving an interstitial vortex will 
sit in the center of the rectangle where the forces from 
the vortices at the pinning sites cancel. For decreasing 
Ly/Lx, the interstitial vortices will have to move through 
a narrower channel between vortices at the pinning sites 
in the y-direction so the interaction between the inter- 
stitial vortices and the pinned vortices will increase. For 
driving along the y direction, for Ly/L^ decreasing from 
1, the critical force is reduced. This can be understood by 
considering that as Lx is reduced the interstitial vortex 
interactions with the pinned vortices in front of the inter- 
stitial vortices and those behind will increasingly cancel. 
As Ly/Lx is increased from 1 the opposite occurs with 
the interaction of the vortices that push on the interstitial 
vortices becoming reduced. 

Experimentally the behavior of the critical currents as 
predicted from equations (5) and (6) can also be used 
to determine whether interstitial vortices are present at 
B/B^ = 2 as opposed to multiple vortices at individ- 
ual pinning sites. In the case of multiple vortices at the 
pinning sites the depinning force would be independent 
of the ratio of Ly/ since the depinning force will not 
be determined by the interactions between the vortices 
due to the symmetry of the overall multi-vortex lattice, 
but will instead be determined from the strength of the 
pinning site only. We also show in Fig. 1 the simula- 
tion results for the square {Ly/Lx — 1.0) and rectangular 
cases [Ly/Lx = 2.0 and Ly/Lx = 0.5) for driving in the 
X and y directions showing excellent agreement with the 
predicted values. 

Another experimental signature of interstitial vortices 
is the strong difference in the critical force between the x 
and y directions at Ly — 2Lx. The critical force in the y 
direction is already more than 50 times higher than the 
critical force in the y-direction. Even experiments with 
a relatively small anisotropy ratio such as Ly/Lx = 1.25 
should see f^/fy « 3. 

These results also show that the maximum pinning can 
be achieved not with square arrays but with rectangular 
arrays. This enhancement is only for one of the directions 
of driving. Figure 1 shows that the depinning force for 
Ly = 2Lx in the y-direction is about 2.5 times higher 
than for the square case. In practice the ratio of Ly/Lx 
has a finite range in which it can be varied due to the 
finite size of the pinning sites. 

We note that wc could calculate analytically the de- 
pinning force only at _B = 2_B0, or for any other B values 
where the interstitial-interstitial vortex interactions effec- 
tively cancel. However, away from the B values leading 



to high symmetry interstitial configurations, the depin- 
ning forces should still be anisotropic. The depinning 
force for an interstitial vortex near a single interstitial 
vacancy will be reduced from its value at B = 2B^ by a 
quantity oc 1 / Lx if driven in the x direction and oc 1 / Ly 
if driven in the y direction since the vortex-vortex force 
goes like oc 1/rij. Similar arguments for an anisotropic 
depinning force can be made for additional interstitial 
vortices added to the B — 2B^ vortex configurations. 

IV. DEPENDENCE OF DEPINNING FORCE ON 
B/B4, IN A SQUARE ARRAY 

In Fig. 2 we show the dependence of the critical de- 
pinning force as a function of B/B^ for a system with 
a square pinning array. Here sharp peaks in can be 
seen at B/ B^ — 1, 2, 3, 4, and 5. Smaller peaks are seen 
for B/Bij, — 6 and 8. At B/B^ — 7 there is no evi- 
dence for a peak. Another interesting feature is that for 
6 < -B / B^ < 8 the critical force remains at an intermedi- 
ate value which is higher than the lowest critical current 
values for B < 5B^. This is a similar trend to that ob- 
served by Metlushko et at who claim that interstitial 
vortices are present for B/B^ > 1 due to the small size of 
the pinning sites. The peaks observed in Fig. 2 are much 
stronger then those observed in magnetization measure- 
ments in simulations with flux-gradient driven vortices 
with short range interactions psf . In that work commen- 
suration enhancements were not seen at B / B^ = 3, 6, or 
7. The vortex configurations observed at the matching 
fields are the same as those for Ref [26]. Another feature 
in Fig. 2 is that the height of the peaks for B < ^B^ 
shows variations with the largest peak at B = AB^ when 
the vortices form a triangular lattice as seen in simula- 
tions and experiments. Some experiments with square 
and magnetic dots have observed strong matching effects 
at every matching field suggesting the presence of 
multi-vortex states in these systems. The inset of Fig. 2 
shows the depinning curve for a system with the same pa- 
rameters as Fig. 2, but with the pinning sites in a random 
arrangement. Here the depinning force decreases with 
increasing field and there are no peaks at the matching 
fields. 

In Fig. 2 commensurability peaks at the fractional 
fields 1/2,3/2, and 5/2 can also be observed. Peaks in 
the critical current at fractional B / B^ have been seen in 
experiments with particularly pronounced peaks at 1/2 
and 3/2 while weaker peaks were seen at 1/4, 1/5, and 
1/16 1^ . The vortex configurations and dynamics at frac- 
tional matching fields are studied in detail elsewhere. We 
do not observe any particular fractional matching except 
at B / B^ = 8.5 where a clear peak in the critical current 
is observed. Such a commensuration peak has not been 
seen in previous simulations or experiments. In section 
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V (c) we show that a stable non-mixing flow state occurs 
at this field. 

V. DYNAMICAL FLOW STATES FOR SQUARE 
PINNING ARRAYS 

A. ID Flow States at B/B^ = 2,4 and 9 

Figure 3(a) shows the flow states just above depinning 
for B /Bij, = 2 where the motion consists of the straight 
ID flow of interstitial vortices between the pinned vor- 
tices. This type of flow was also seen in simulations at the 
same field for vortices with bulk interactions. In Fig. 3(b) 
for B/B^ = 4 a similar ID straight vortex flow is ob- 
served. Here the vortices at the pinning sites also remain 
pinned and there are additional interstitial vortices lo- 
cated between the pinning sites in the x direction that 
also remain pinned. The mobile interstitial vortices can 
move in a straight unhindered ID path while the immo- 
bile interstitial vortices cannot move in ID paths without 
entering a pinning site or coming close to the vortices lo- 
cated in the pinning sites. In Fig. 3(c) at B/B^ — 9 
the same type of ID interstitial flow as in Fig. 3(a) and 
Fig. 3(b) is observed but in this case there are two mobile 
rows of vortices between the pinning sites and two immo- 
bile interstitial vortices between vortices at the pinning 
sites. In Fig. 3(a,b,c) the vortex motion can be seen to be 
elastic with respect to the mobile vortices in which the 
moving vortices keep their same moving neighbors. Fur- 
ther, the vortices always flow in the same paths. At fields 
with a fraction of 0.05 higher or lower than the matching 
fields in Figs. 3(a,b,c) the initial vortex motion occurs at 
the location of the extra vortex or vacancy in the ordered 
interstitial lattice. The depinning occurs at a lower driv- 
ing force than that at which the commensurate vortex 
configuration depins. The flow of these extra interstitials 
or vacancies will again be in a ID path along the direc- 
tion of drive; however, the motion is not continuous but 
occurs in a soliton or pulse fashion with the extra inter- 
stitial or vacancy exchanging places with pinned vortices 
as it propagates. A similar soliton like motion of vortices 
along the pinning sites has also been observed previously 
in simulations p3| and experiments [^. 

B. Sinusoidal Flow States At B/B4, = 5,8, and 8.5 

In Fig. 4(a) we show the flow states at B/ B^ — 5 which 
shows sinusoidal flow states of the interstitial vortices. 
Here the vortex motion is not strictly ID in the direction 
of drive but shows a periodic motion in the transverse 
direction as well. In addition the vortices flow in the 
same paths and vortices from one channel do not mix 
with vortices in the other channels. In Fig. 4(b) a simi- 
lar sinusoidal flow as seen at B/B^ = 5 is observed with 



an additional square sub-lattice of pinned interstitial vor- 
tices at the center of the pinning plaquette similar to the 
vortex configuration at B/B^ ~ 2. For higher drives the 
immobile interstitial vortices depin and the vortex lattice 
enters a new flow phase. In Fig. 4(c) for B/B^ = 8.5, 
where a small peak in the critical current is seen, an in- 
teresting flow state is observed where there is a pinned 
interstitial lattice for every other interstitial vortex in 
the middle of the pinning array plaquette. Additionally 
there are two immobile vortices between pinning sites in 
the x-direction. The vortex flow is a combination of ID 
flow and sinusoidal flow with the sinusoidal flow occurring 
at every other plaquette containing a pinned interstitial 
vortex. 

C. Bistable Flow States in the Vicinity of B/B^ = 3 

Fig. 5(a,b) shows typical vortex flow patterns observed 
for fields 2.75 < B/B^ < 3.5. Here the interstitial 
vortices move in periodic meandering paths around the 
pinned vortices. Interestingly the direction of the vor- 
tex flow is not in the direction of the driving force but 
at 45 degrees in either transverse direction. Since either 
direction is equivalent the flow jumps between the two 
states as seen in the measured transverse Vy velocities. 
In Fig. 6, Vy shows a small amplitude periodic component 
caused by the winding nature of the vortex flow, along 
with low frequency large amplitude jumps of Vy from 
positive to negative indicating that the net vortex mo- 
tion direction is changing. For different system sizes the 
characteristic of the flow paths remains the same. Exper- 
imentally this flow state can be observed with transverse 
noise measurements. For stronger drives there is a tran- 
sition to a more disordered flow in the direction of the 
driving force. 

D. Disordered Vortex Flow States and Noise 

In Fig. 7(a,b,c) we show the vortex flow states at var- 
ious incommensurate fields showing varying degrees of 
disordered or mixing flow states. In Fig. 7(b) the vor- 
tex trajectories for B/B^ = 4.95 are plotted showing 
that the general features of the stable channel flow from 
B/Bfj, = 5 are still present. There are, however, now 
some vortices that can be seen to jump between adjacent 
channels. In Fig. 7(a) for B/B^ = 4.5 the vortex trajec- 
tories are more disordered and no discernible preferred 
channels are visible. This same type of flow is seen for 
B/B^ = 5.5. There are regions near the pinning sites 
where repulsion from the vortices located in the pinning 
sites keeps other vortices from approaching closer than 
a certain distance. In Fig. 7(c) we show the vortex tra- 
jectories for B/B^ = 7 showing that the disordered flow 
is again similar to what was seen at B/B^ = 4.5. We 
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note that this is the matching field at which a peak was 
absent in the depinning force. In general we observe that 
the non-mixing flows occur only near certain integer or 
fractional matching fields where the initial pinned vortex 
lattice has an ordered state. The flows then become in- 
creasingly disordered as the fields move away from these 
matching configurations. 

To further compare the experimentally measurable 
noise signals of the ordered and disordered vortex states 
we plot in Fig. 8 the time series and the correspond- 
ing Fourier transforms (FT's) of the average vortex ve- 
locity Vx for a fixed applied driving force. Each data 
set is analyzed for 30,000 MD steps. In the figures we 
only show a portion of the time series for clarity. For 
B/B^ = A (Fig. 8(a,b)), where the flow was seen to occur 
in ID channels as seen in Fig. 3(b), only a single peri- 
odic component and higher harmonics can be seen. For 
B/B^ — 4.5 (Fig. 8(c,d)) where the flow was disordered 
as seen in Fig. 7(a), a more random noise signal is ob- 
served and the FT shows a broad spectra. At B /B^ = 5 
(Fig. 8(e,f)) where an ordered flow state was observed as 
seen in Fig. 4(a), a clear single frequency signal is again 
observed. In Fig. 9 we show the time series and FT for 
B/B^ — 3 where more complicated vortex trajectories 
were observed as seen in Fig. 5(a,b). Here the signal 
shows several different periodicities. In general we find 
periodic signals where orderly or partially ordered flows 
are observed, and broad noise spectra where incommen- 
surate or disorderly flows occur. 

E. Longitudinal and Transverse Noise 

In Fig. 10(a,b) we show the time series for both the 
sum of the velocities in the a;-direction and j/-direction 
simultaneously for (a) B/B^ = 4 and (b) B/Bff, = 5. For 
B /B^ = 4 longitudinal velocities show the periodic signal 
as shown earlier while the transverse velocities are zero 
as the vortices are moving strictly in the x-direction in 
ID channels as seen in Fig. 3(b). For B/B^ = 5 both 
the longitudinal and transverse velocities show a pro- 
nounced periodic signal as the vortex trajectories show 
a sinusoidal flow as seen in Fig. 4(a). These results sug- 
gest that the features of the flow, such as whether it is 
winding in the transverse direction or strictly ID, can 
be probed experimentally with transverse noise measure- 
ments. The periodic time signals at the matching fields 
can also be probed with an applied AC drive superim- 
posed on a DC drive. Shapiro steps can be observed, as 
recently determined in experiments and simulations at 
B/Bfj, = 2. Where the vortex flows are highly ordered 
one would then expect Shapiro steps |2^ ], while at the 
disordered flow states the phase-locked steps would be 
absent or strongly reduced. 



VI. DEPENDENCE OF DEPINNING FORCE ON 
FIELD FOR RECTANGULAR PINNING ARRAYS 

A. Rectangular Pinning Array with Lx/Ly — 2 

In Fig. 11 (a, b) we show the dependence of the critical 
depinning force on B/B^ for systems with similar pa- 
rameters as studied for the square array with Lx = 2Ly 
for driving in the x-direction and y-direction. Here is 
much higher than fy for all B > . For B / B^ < 6 com- 
mensurability peaks are observed for driving in both the 
X and y directions; however, the peaks in the z-direction 
do not show the sharpness of the peaks seen for the y- 
direction. The peaks for both driving directions are much 
broader for B/B^ > 5. There are no clear peaks at 
B/B^ = 7 or 9; however, there is some evidence for a 
drop-off in the depinning force right after these fields. 

B. Vortex Configurations for Lx/ Ly = 2 

In Fig. 12(a-i) we show the vortex configurations for 
the various matching fields. The vortex configurations at 
B/Bij, = 1 and 2 (Figs. 12(a,b)) are rectangular in or- 
der. For B/B^ = 3 (Fig. 12(c)) the vortex crystal has a 
diagonal ordering. At B/B^ ~ A (Fig. 12(d)) a square 
vortex lattice is stabihzed. At B/B^ ~b (Fig. 12(e)) the 
vortex lattice does not form a simple rectangular lattice 
but an ordering is still present as indicated by the unit 
cell. We note that at B/B^ = 5 a peak or anomaly in 
the depinning force is not observed. For B/B^ = 6 and 

8, ordered distorted hexagonal vortex lattices are stabi- 
lized. At B/B^ = 7 and 9 the vortex lattice is disordered 
which again corresponds to there being no peaks in the 
critical depinning force. These results indicate that at 
the matching fields where an ordered lattice with a sim- 
ple rectangular or hexagonal ordering can occur, an en- 
hancement in the critical current can be observed. 

VII. DYNAMICAL FLOW STATES FOR 
RECTANGULAR PINNING ARRAYS 

In Fig. 13(a-f) we show the flow states at depinning 
in the rectangular arrays for B/B^ = 2,3,4,5,6, and 

9, respectively, for driving along the long side or x- 
direction. ID flows that are similar to those seen in the 
square arrays occur for B/B^ — 2 and 6 (Fig. 13(a) and 
Fig. 13(e)). Disordered flow phases occur for B/B^ — 3, 
4 and 9 (Fig. 13(b), Fig. 13(c) and Fig. 13(f)) where 
only minimal anomalies were seen in the depinning force. 
A disordered phase is also observed for B/B^ = 7 (not 
shown) where there was again no anomaly in the de- 
pinning force. An ordered flow phase is observed for 
B/Brj, = 8 which is described in the next section where 
a peak in the depinning is observed. At B/B^ = 5 
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(Fig. 13(d)). where a static ordered lattice was stabilized 
but an anomaly in the depinning force was not observed, 
we find an interesting periodic winding flow of interstitial 
vortices. Also pairs of immobile interstitial vortices are 
located between every other pair of pinning sites. 

A. Pinched Sinusoidal Flow For B/B4, > 5 for Driving 
Along the Short Edge 

A new type of stable flow phase that appears for the 
rectangular arrays for B/B^ > 4 at higher drives is a 
pinched sinusoidal flow phase such as seen in Fig. 14(d) 
for B = 8. Wc follow the development of this phase in 
Fig. 14 for B/B0 = 8 where the initial depinning occurs 
through the almost ID flow of interstitial vortices where 
two interstitial vortices remain pinned between the pin- 
ning sites as seen in Fig. 14(a) where fd = 0.1. There is 
additional slower periodic wandering of the moving vor- 
tices in the y-direction that broadens the channel. This 
periodic transverse wandering of the channels can also 
be seen to push the immobile interstitial vortices in the 
transverse direction. The trajectories from the longitu- 
dinally moving vortices are pinched where the moving 
vortices are directly adjacent to the pinning sites. In 
Fig. 14(b) with fd = 0.2 we show that for increased driv- 
ing forces there is a transition to a disordered state as 
the immobile interstitial vortices begin to depin. For 
fd = 0.3 (Fig. 14(c)) the flow starts to reorganize to the 
pinched sinusoidal flow which becomes fully developed 
for fd = 0.325 as seen in Fig. 14(d). 

In the pinched sinusoidal flow the vortices show a sta- 
ble channel flow similar to the sinusoidal flows seen in 
the square arrays; however there is a pinched feature as 
the vortices flow adjacent to the pinning sites. As the 
vortices move through the pinched area only one vortex 
goes through the intersections at a time with a vortex 
from the upper and lower part of the branches alternat- 
ing. For the other matching fields for B/B^ > 4 the 
same type of flow patterns are observed. The same flow 
patterns are observed at fields not too far from commen- 
suration with the mobile rows carrying a varying number 
of vortices. 

B. Additional Dynamical Flow Phases for B/Btj, < 4 

In Fig. 15(a) we show an interesting type of dynami- 
cally induced rotational motion ordered flow phase that 
occurs at and near B /B^ = 2.5 for driving along the long 
side of the pinning lattice. The depinning is due to the ID 
flow of a portion of the interstitial vortices which move 
between every other pinning row. Another portion of in- 
terstitial vortices cannot pass through the pinned vortices 
and show a remarkable paired collective rotational mo- 
tion. The rotations occurs as a pair of interstitial vortices 



rotate in opposite directions as the ID moving intersti- 
tial vortices pass. Additionally every other rotating pair 
is out of phase with the adjacent rotating pairs. The ro- 
tational motion occurs due to the fact that the rotating 
interstitial vortices are located in a shallow potential well 
created by the vortices located at the pinning sites. As 
the ordered ID moving vortices move past they push the 
interstitials vortices that are not moving longitudinally. 
Although these pairs do not show a net motion in the x- 
direction they are still taking part in the dissipation and 
slow down the ID moving vortices. For higher drives the 
rotating interstitials can be dislodged firom their orbits 
and the flow becomes more disordered. 

In Fig. 15(b) we show the ordered flow phase that oc- 
curs for B/B^ = 4.5 where the initial depinning occurs 
through the sinusoidal flow of a portion of the interstitial 
vortices. Unlike the previous sinusoidal flows observed 
in this case the channels weave across the pinning sites 
rather then strictly between two rows of pinning sites. 
For higher drives the other immobile vortices depin and 
the flow enters a disordered phase. 

C. Flow Phases For Driving Along Long Edge 

In Fig. 16(a-f) we show the flow states just above de- 
pinning for driving along the short side of the rectangular 
pinning array for B /B^ = 2, 3, 4, 7, 8, and 9, respectively. 
The motion a,t B = 2B^ (Fig. 16(a)) is ID where a sin- 
gle channel of interstitials is moving while for B /B^ = 8 
(Fig. 16(e)) three rows of vortices are flowing and an im- 
mobile interstitial vortex is located between the pinning 
sites. An almost ID flow is seen at B/B^ = 7 (Fig. 16(d)) 
in which 3 rows of interstitial vortices are moving. We 
note that for B / B^ = 7 a peak or matching anomaly was 
not present. The vortex configurations for B / B^ = 7 
shows that along the immobile row of vortices there are 
16 vortices witli 8 located in the pinning sites while the 
mobile rows each contain 13 vortices so that the mov- 
ing vortices are not commensurate with the periodicity 
of the potential created by the immobile vortices. For 
the B/Bfj, = 8 matching field where a peak was seen, 
the vortex configurations again show 16 vortices along 
the immobile row but the moving rows each contain 16 
vortices which can be commensurate with the immobile 
vortices. The flow paths for B/B^ = 3 (Fig. 16(b)) and 
B/B(j, — 4 (Fig. 16(c)) each show two rows of moving 
interstitial vortices in a sinusoidal flow with the flow at 
B/B^ = 4 showing a a larger amplitude. For B/B^ = 9 
(Fig. 16(f)) where there was an absence of the peak in 
the depinning force the flow is mostly disordered with 
some remnant of the flow phase seen at B/B^ — 8. 

In Fig. 17(a,b) we show the flow phase at B/B^ = 
5, and 6, respectively, for driving along the short edge. 
The flow at B / B^ = 7 shows a remarkable braiding flow 
where the moving interstitial vortices flow in a crossing 
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pattern. A moving vortex that starts in the region almost 
underneath the pinning sites wiU move at an angle until it 
reaches the opposite side of the channel two pinning sites 
up while other vortices cross in the opposite direction. 
At B/ = 6 a similar flow to that seen a.t B/ B^ = 5 is 
observed; however, the motion is not perfectly ordered. 

VIII. CONCLUSION 

We have analyzed numerically the vortex pinning and 
dynamics in square and rectangular pinning arrays for 
thin film superconductors. We predict an anisotropic 
critical current for the rectangular pinning arrays and 
calculate analytically the ratio of the critical depinning 
forces for different aspect ratios at B/B^f, = 2. These 
results show that as the critical current is enhanced in 
one direction it is reduced in the other direction. The 
maximum critical current can be achieved for driving in 
certain directions of the rectangular pinning array. Nu- 
merical simulations also find that the anisotropy in the 
critical depinning force occurs for B/B^ < 4 and is less 
pronounced for higher fields. In the case of the square 
pinning arrays we observe pronounced commensurabil- 
ity effects at most integer matching fields and a missing 
matching peak at B/B^ = 7 in good agreement with 
recent experiments. 

In the rectangular pinning arrays with an aspect ra- 
tio of 2:1 we observe an anisotropic depinning threshold 
in which the easy flow direction shows an overall lower 
depinning force but sharper and more pronounced com- 
mensuration effects while in the other direction the over- 
all depinning force at both the commensurate and non- 
commensurate fields is higher, but the matching effects 
are considerably reduced. We observe integer matching 
effects in both driving directions up to B/B^ = 5, after 
which only every other matching field shows an enhanced 
critical current. The vortex configurations show that at 
the matching fields where enhanced critical currents are 
observed an ordered vortex arrangement is formed while 
at the matching fields where there is no enhancement the 
vortex arrangement is disordered. 

For both the square and rectangular arrays we find that 
a remarkable variety of intricate dynamical flow phases 
can be realized and in general two classes of flow phases 
can be observed: stable channel flow, where vortices flow 
in the same paths in identifiable channels and vortices 
from one channel do not mix with vortices in another 
channel; and a disordered or mixing flow, where vortices 
from different channels mix or no identifiable channels 
occur. The particular fiow states we observe include a 
ID fiow of interstitial vortices. For the square arrays at 
B/B^ = 5,6, 8 and 8.5 a sinusoidal flow occurs where a 
portion of the interstitial vortices are moving in winding 
paths while another portion of interstitial vortices remain 
pinned along with the vortices at the pinning sites. For 



increasing drives the immobile interstitial vortices can 
depin and the vortices can enter a more disordered fiow 
phase. Near B/B^ = 3 a. remarkable bistable flow phase 
is observed where the vortex flow is not in the direction 
of drive, but is alternating from the +45 degree direc- 
tion to the —45 degree direction. This motion can also 
be seen in the transverse vortex velocities. Away from 
commensurability these ordered fiow phases become in- 
creasingly disordered. We also show that these phases 
can be probed experimentally with noise spectra. In or- 
dered phases the noise spectra shows distinctive narrow- 
band features, while the disordered flows show a broad 
spectra. In addition we show that with transverse noise 
measurements the ID and sinusoidal flows can be dis- 
tinguished where the sinusoidal flows will show a narrow 
band transverse signal. 

In the rectangular arrays several of the phases observed 
in the square arrays are also observed. A new ordered 
flow phase that occurs for driving in the short direction 
for a wide range of fields for high driving is a braided 
channel flow. At the matching fields the braided flow 
is uniform with each channel carrying the same number 
of vortices while at incommensurate fields the channels 
can carry different numbers of vortices. The signature of 
the onset of this braided flow can be seen in the current- 
voltage characteristics. 

Acknowledgments 

We thank C.J. Olson for critical reading of this 
manuscript, and I.K. SchuUer, J. Martin, L.E. DeLong, 
L. Van Look, V. Metlushko, V.V. Moshchalkov, A. Hoff- 
man, F. Nori, and R.T. Scalettar for useful discussions. 
This work was supported by the Director, Office of Ad- 
vanced Scientific Computing Research, Division of Math- 
ematical, Information and Computational Sciences of 
the U.S. Department of Energy under contract number 
DE-AC03-76SF00098 as well as CLC and CULAR (Los 
Alamos National Laboratory). 



[1] A.T. Fiory, A.F. Hebard, and S. Somekh, 

Appl. Phys. Lett. 32, 73 (1978). 
[2] V. Metlushko et ai, Solid State Commun. 91, 331 (1994). 
[3] M. Baert et ai, Phys. Rev. Lett. 74, 3269 (1995); 

V.V. Moshchalkov et al, Phys. Rev. B 54, 7385 (1996); 

V.V. Moshchalkov et al, Phys. Rev. B, 3615 (1998); 

V. Metlushko et al, Phys. Rev. B 59, 603 (1999). 
[4] M. Baert et al, Europhys. Lett. 29, 157 (1995). 
[5] J.Y. Lin el al, Phys. Rev. B 54, R12714 (1996); 

A. Castellanos et al, Appl. Phys. Lett. 71, 962 (1997). 
[6] R. Surdeanu, R.J. Wijngaarden, R. Wordenweber, and 

R. Griessen, to be published. 



8 



[7] 

[8] 

[9] 
[10] 

[11] 
[12] 

[13] 



[14] 

[15] 
[16] 



[17] 

[18] 
[19] 
[20] 
[21] 
[22] 
[23] 

[24] 

[25] 

[26] 

[27] 

[28] 

[29] 

[30] 

[31] 



[32] 

[33] 
[34] 

[35] 
[36] 

[37] 

[38] 

[39] 



and B. Pannetier, 



A. Bezryadin, Yu. B. Ovchinnikov, 
Phys. Rev. B 53, 8553 (1996). 
E. Rosseel et al, Phys. Rev. B 53, R2983 (1996). 
K. Harada et al , Science 271, 1393 (1996). 
S. Field 



ond-mat/0003415; A.N. Grigorenko et 



et al.. 
al, to be published. 

L. Van Look et al, Phys. Rev. B 60, R6998 (1999). 
V. Metlushko et al, Phys. Rev. B 60, R12 585 (1999). 
J.I. Martm et al., Phys. Rev. Lett. 79, 1929 (1997); 
Y. Jaccard et al., Phys. Rev. B 58, 8232 (1998); A. Hoff- 
man, P. Prieto and I.K. SchuUer, Phys. Rev. B 61, 6958 
(2000). 

D.J. Morgan and J.B. Ketterson, Phys. Rev. Lett. 80, 
3614 (1998). 

Y. Fasano et al, Phys. Rev. B 60, R15047 (1999). 
A. Terentiev, D.B. Watkins, L.E. De Long, D.J. Morgan 
and J.B. Ketterson, Physica C 324, 1 (1999); A. Teren- 
tiev et al, Phys. Rev. B 61, R9249 (2000). 
V.V. Metlushko et al., Europhysics Letters 41, 333 
(1998). 

J.L Martm et al, Phys. Rev. Lett. 83, 1022 (1999). 

J.L Martm et al, Phys. Rev. B 62, 9110 (2000). 

L. Van Look et al, to be published. 

S. Kolesnik et al., to be published. 

V.V. Metlushko et al., to be published. 

C. Reichhardt, C.J. Olson, and F. Nori, Phys. Rev. B 57, 

7937 (1998). 

C. Reichhardt, C.J. Olson, and F. Nori, 
Phys. Rev. Lett. 78, 2648 (1997). 

C. Reichhardt, R.T. Scalettar, G.T. Zimanyi, and 
N. Gr0nbech- Jensen, Phys. Rev. B Rll 914 (2000). 
V.I. Marconi and D. Dommgue z, Phys. Rev. Lett. 82, 



Phys. Rev. Lett. 70, 
Bhattacharya, Phys- 
Hellerqvist et al. , 
W. Henderson et al. 



4922 (1999); |cond-mat/000812E 
K.D. Fisher, D. Stroud, and L. Janin, Phys. Rev. B 60, 
15371 (1999). 

A. Pruyjmboom et al, Phys. Rev. Lett. 60, 1430 (1988); 

M.H. Theunissen et al, ihtd, 77, 159 (1996). 

S. Anders et al, Physica C 332, 35 (2000); and to be 

published. 

R. Besseling, R. Niggebrugge, and P.H. Kes, 
Phys. Rev. Lett. 82, 3144 (1999). 
S. Bhattacharya and M.J. Higgins, 
2617 (1993); M.J. Higgins and S. 
ica C 257, 232 (1996); M.C. 
Phys. Rev. Lett. 76, 4022 (1996); 
Phys. Rev. Lett. 77, 2077 (1996). 

A.C. Marley, M.J. Higgins and S. Bhattacharya, 

Phys. Rev. Lett. 74, 3029 (1995). 

U. Yaron et al. Nature 376, 753 (1995). 

A. Duarte et al, Phys. Rev. B 53, 11336 (1996); F. Pardo 

et al, Phys. Rev. Lett. 78, 2644 (1997). 

F. Pardo et al, Nature 396, 348 (1998). 

A.M. Troyanovski, J. Aarts, and P.H. Kes, Nature 399, 
665 (1999). 

A.E. Koshelev and V.M. Vinokur, Phys. Rev. Lett. 73, 
3580 (1994). 

K. Moon, R.T. Scalettar, and 

G. T. Zimanyi, Phys. Rev. Lett. 77, 2278 (1996); S. Ryu 
et al, Phys. Rev. Lett. 77, 5114 (1996). 

H. J. Jensen, A. Brass and 
A.J. Berlinsky, Phys. Rev. Lett. 60, 1676 (1988); A.Ch- 



Shi and A.J. Berlinsky, Phys. Rev. Lett. 67, 1926 (1991); 
N. Gr0nbech- Jensen, A.R. Bishop, and D. Dommguez, 
Phys. Rev. Lett. 76, 2985 (1996); M.C. Faleski, 
M.C. Marchetti, and A.A. Middleton, Phys. Rev. B 54 
12427 (1996); S. Spencer and H.J. Jensen, Phys. Rev. B 
55, 8473 (1997); D. Dommguez, Phys. Rev. Lett. 82, 
181 (1999); A.B. Kohon, D. Domi'nguez, N. Gr0nbech- 
Jensen, Phys. Rev. Lett. 83, 3061 (1999). 
[40] C.J. Olson, 
C. Reichhardt, and F. Nori, Phys. Rev. Lett. 81, 3757 
(1998). 

[41] T. Giamarchi and P. Le Doussal, Phys. Rev. Lett. 76, 
3580 (1996); 78, 752 (1997); P. Le Doussal and 
T. Giamarchi, Phys. Rev. B 57, 11 356 (1998); 
L. Balents, M.C. Marchetti, and L. Radzihovsky, 
Phys. Rev. Lett. 78, 751 (1997); Phys. Rev. B 57, 7705 

(1998) ; S. Scheidl and V.M. Vinokur, Phys. Rev. E 57, 
2574 (1998). 

[42] N. Gr0nbech- Jensen, Int. J. Mod. Phys. C 7, 873 (1996); 
N. Gr0nbech- Jensen, Comp. Phys. Comm. 119, 115 

(1999) . 



0.10 



0,01 \L 





fdiiy ^ — — - 




fjlx \. : 


-/ 




/ , , , , 





0.5 



l.( 



2.0 



2.5 



1.5 

FIG. 1. Critical depinning force for driving in the 
a;-direction and y-direction vs the aspect ratio of the rect- 
angular pinning array for B — 2B^ as predicted by Eq. (5) 
and Eq. (6) for a constant Ly — 2.0 with Lx varied from 4.0 
to 1.25. The filled circles are the results from simulations. 
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FIG. 2. Critical depinning force vs the field B/B^ in a 
system with a square pinning array for nLx = nLy = 12.0. 
Peaks in the depinning force can be seen at most of the match- 
ing fields with a clear missing peak at B/B^ — 7. In addition 
some clear peaks can be seen at the fractional matching fields 
B/Bijy = 1/2, 3/2 and 5/2. The inset shows the depinning 
force vs B/B^ for a system with the same parameters but 
with the pinning sites in a random arrangement showing the 
absence of matching effects. 
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FIG. 3. Vortex trajectories for the system in Fig. 2 just 
above depinning for (a) B/B^ = 2, showing the ID flow of 
single rows of interstitial vortices, (b) B/B^ — 4, showing the 
ID flow of 2 rows of vortices between pinning sites. Between 
each pair of pinning sites is a single immobile interstitial vor- 
tex, (c) B/Brj, = 9, showing the ID flow where 3 rows of 
vortices flow between the pinning sites while a pair of immo- 
bile interstitial vortices are located between the pinning sites. 
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FIG. 4. Vortex trajectories just above depinning for (a) 
B/B^ = 5 showing that the flow occurs with two periodic 
winding rows of mobile vortices moving between each pinning 
row. {h)B/B^ = 8 showing a more sinusoidal like flow where 
again two rows of vortices move between the pinning rows. 
In addition, in the center of each pinning plaquette is an im- 
mobile interstitial vortex, (c) B/B^ = 8.5 shows a pinched 
flow. An immobile interstitial vortex is located in the center 
of every other plaquette and a pair of immobile interstitial 
vortices are located between every other pair of pinning sites. 
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FIG. 6. Transverse, or y-direction, vortex velocities for the 
system in Fig. 5 where the driving is in the x-direction. The 
vortex motion locks into the ±45° direction for 150000 MD 
steps before switching to the other direction. The additional 

smaller scale periodic component is due to the winding nature 
of the vortex flow in the channels as seen in Fig. 5. 





FIG. 5. (a) Vortex trajectories for B/B^ = 3 showing 
tli(^ periodic winding vortex motion at +45°. (b) The vortex 
trajectories for the same system and time interval as in (a) 
after the vortex motion switches to the —45°. 
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FIG. 7. Vortex trajectories for (a) B/B^ = 4.5 showing 
a disordered flow where the trajectories change over time. 
Vortices do not move through the regions near the occu- 
pied pinning sites due to the vortex-vortex repulsion, (b) 
B/Bcp — 4.95, showing that flow is more ordered with the 
same features of the periodic winding channels observed at 
B/B^ — 5. If the trajectories are drawn for a longer time 
the plots become increasingly smeared and the plot will be- 
come indistinguishable from that seen in (a), (c) The vortex 
trajectories for the same time interval as in (a) and (b) for 
B/B^ — 7 showing a disordered flow pattern. 
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FIG. 8. Time series of the sum of longitudinal vortex 
velocities and Fourier transforms for the vortex flow. (a,b) 
B/B^ = 4 where a periodic signal occurs due to the ID flow 
of vortices (Fig. 3(b)) through the periodic potential created 
by the pinned vortices. (c,d) B/B^ = 4.5 where a broad 
spectra occurs due to the disordered flow as seen in Fig. 7(a). 
(e,f) B/B^ — 5 where again a periodic signal occurs as the 
vortices move in periodic winding channels as seen in Fig 4(a). 
The FT's were taken on data sets of 30000 MD steps. 
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FIG. 9. Time series of the sum of the longitudinal vor- 
tex velocities and Fourier transform for the vortex flow at 
B/B^ — 3, where the flow occurs in winding channels. Here 
several different periodicities can be distinguished. 
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FIG. 10. The time scries for the sum of the longitudinal 
velocities 14 and transverse velocities Vy for (a) B/B^ = 4 
and (b) B/B^ = 5 for a fixed drive of fa/ fp = 0.4. In (a) 
where the vortex motion occurs by the ID flow of interstitial 
vortices as seen in Fig. 3(b) the longitudinal velocities show a 
periodic component while the transverse velocity is zero. For 
(b) where the vortex motion was sinusoidal both the x and y 
components of the velocities show a periodic signal. 
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FIG. 11. Critical depinning forces versus B/B^ for driving 
in the x and y directions for a system with a rectangular pin- 
ning array with L^/Ly — 2. (a) Dashed curve is the depinning 
line for driving in the x-direction and the solid curve is the 
depinning line for driving in the y-direction. For B/B^ < 4 
there is a clear anisotropy in the depinning force. The driv- 
ing along the x-direction gives a higher depinning force ex- 
cept at B/B^ = 1, where the depinning forces both equal 
fp. In (b) the same curve is plotted with the x-direction 
driving curve (dashed-line and open circles) shifted up for 
clarity. The matching peaks for the a;-diroction driving are 
strongly reduced as compared to the depinning curve for driv- 
ing in the y-diroction. Missing peaks for both curves occur at 
B/B^ = 5, 7, and 9. 
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FIG. 12. Vortex positions (black circles) and pinning 

sites (open circles) after annealing from a higher tempera- 
ture state in a rectangular pinning array with L^/Ly =2 for 
(a) B/B^ = 1, (b) 2, (c) 3, (d) 4, (e) 5, (f) 6, (g) 7, (h) 8, 
and (i) 9. 
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FIG. 14. Vortex trajectories for driving in the a;-direction 
in a rectangular pinning array with L^/Ly =2 at B/B^ = 8. 

(a) The flow patterns for fd/ fp = 0.2, where the flow is or- 
dered and a portion of the interstitial vortices are immobile. 

(b) fd/ fp = 0.3 flow becomes disordered, (c) fd/ fp = 0.35 
flow begins to reorganize, (d) fa/ fp = 0.37 flow is ordered 
with all the interstitial vortices moving. 
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(a) (b) 
FIG. 15. Vortex trajectories for driving in the a;-direction 
in a rectangular pinning array with L^/Ly = 2 just above 

depinning. (a) B/B^ — 2.5; Vortex motion consists of the ID 
flow of interstitial vortices between every other pinning row 
and the rotational motion of pairs of interstitial vortices, (b) 
B/B^ = 4.5; Motion consists of winding channels that weave 
between the pinning sites. 
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(e) (f) 

FIG. 16. Vortex trajectories for driving in the y-direction 
in a rectangular pinning array with L^/Ly = 2, just above 
depinning. (a) B/B^ = 2, (b) 3, (c) 4, (d) 7, (e) 8 and (f) 9. 




(a) (b) 
FIG. 17. Vortex trajectories for driving in the y-dircctiori 
in a rectangular pinning array with L^/Ly = 2, just above 
depinning. (a) B/B^ = 5 shows an intricate braiding flow 
pattern, (b) B/B^ = 6 shows a similar pattern as in (a) but 
with some disorder. 
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